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1. Introduction 

The study of the asymptotic properties of random matrix models in the limit of large 
matrix size has revealed deep connections between random matrix theory and the 
theory of integrable systems. The origin of these connections is nowadays much better 
understood thanks to the use of Riemann-Hilbert (RH) techniques, which are one of 
the main tools in the inverse scattering theory of integrable systems [TJ [2j [3j HJ IS]- 
Indeed, the occurrence of RH problems leads naturally to Lax pair equations and to the 
corresponding string equations [6J, which in turn explain why integrable equations of 
Painleve type arise in the double scaling limit of critical random matrix models [7]. 
The main goal of this work is to describe in detail and in a general setting (i.e. 
independently of particular examples which will be used only for illustrative purposes) 
the correspondence between the full Painleve I and Painleve II hierarchies and critical 
random matrix models. 

Specifically, we consider random Hermitian matrix models with partition function 



:,v / o - -' ' ' ' Y[( x i - xj) 2 dxi- ■ -dx N , (1) 



where the integration is performed over the eigenvalues x±, . . . , Xjy and V is a polynomial 
potential of even degree and real coefficients (coupling constants) g = (gx, . . . ,g2 P ), 

2p 

V(z) = J29nZ n , g 2p >0. (2) 

n=l 

The large N asymptotics of a matrix model is related to the large N asymptotics of 
the orthogonal polynomials P n ^{x) with respect to the exponential weight e~ NV ( x \ i.e. 
monic polynomials 

Pn, N (x) = x n + a n _ 1 x n - 1 + --- (3) 

that satisfy 

/oo 
P k , N (x)PiAx)e~ NV(x) dx = &kih k , N . (4) 
-oo 

In fact there are several methods [H [9], [TO], [TTJ [12] that reduce the calculation of the 
large iV asymptotic behavior of the free energy 



iV 



-N~ 2 In Z N (5) 



to the calculation of the asymptotic behavior of the recurrence coefficients r n ^ N and s n ^ 
in the three-term recursion relation 

xP ntN (x) = P n+ i,jv(x) + S ntN P ntN (x) + r n)N P n _i, N (x) . (6) 

The structure of these asymptotic expansions depends on the number of disjoint 
intervals (cuts) in the support J of the asymptotic eigenvalue density p(x). In fact, the 
asymptotic expansions in the multi-cut (two or more cuts) case are not in general simple 
power series in iV -1 but have a complicated quasi-periodic dependence on N [4ll5l[T3|IT4]. 
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However, for even potentials the recurrence coefficient s n ^ is identically zero, and if 
in addition the potential is in the regular (noncritical) two-cut case with eigenvalue 
support J = (— j3, —a) U (a, (3), then a substantial simplification of the general behavior 
occurs [13]: the odd terms of the recurrence coefficient r 2n +i,N and the even terms r 2n ,N 
admit (different) asymptotic expansions whose leading terms can be gathered in a single 
formula, 



Our goal is to study the critical behavior of these matrix models with respect 
to a temperature parameter T > 0, i.e. we consider the family of models Zn(T) with 
potentials V(z)/T or, equivalently, with coupling constants g/T. Critical matrix models 
correspond to points T c where the asymptotic free energy ceases to be analytic as a 
function of T. Some of these critical models are related to conformal (p, q) minimal 
models [74 [IS QiE EE] and have asymptotic eigenvalue densities with rational singularities 



near one of the endpoints a of J. It has been known for a long time in the physics 
literature that the double scaling limit of critical matrix models in the one-cut and two- 
cut cases reveals the presence of equations belonging to the Painleve I and Painleve II 
hierarchies respectively (see [18] for a description of these hierarchies). However, 
most of the existing work, specially for the two-cut case, deals only with particular 
examples [5j [7J [121 EE 1201 [ZH [22] and does not offer a general characterization of the 
specific member of the Painleve hierarchy associated to a critical model in terms of the 
coupling constants g and of the critical value of the recurrence coefficient. More recent 
works [T6l [T7] do prove the occurrence of the full Painleve hierarchies in the context of 
critical matrix models, but again they do not give an explicit correspondence between 
Painleve equations and critical models. 

To arrive at this general correspondence, we first perform a detailed analysis of the 
large N expansions of the recurrence coefficient r n ^r for general even potentials in the 
one-cut and two-cut cases, with special emphasis on the latter. Our analysis is based on a 
method for solving continuum limits of the discrete string equation [7J. This method uses 
the resolvent of the Lax operator of the underlying Toda hierarchy and can be applied 
to obtain the large N expansions both in regular and in critical models [231 El] • The 
asymptotic behavior of r n jv f° r regular models is given by power series in N~ 2 whose 
coefficients can be determined recursively. For critical models we apply the double 
scaling method to calculate a family of asymptotic solutions of the string equation 
(symmetric solutions) which are series in fractional powers of N' 1 whose coefficients 
are constrained by ordinary differential equations in the scaling variable. It is at this 
point when we are able to give the precise connection between the complete Painleve I 
and Painleve II hierarchies and certain families of critical models. In particular we show 
how the Painleve II hierarchy is connected to critical models featuring the merging of 




(7) 




(8) 



two cuts [TT1|2T1|22]. 
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The paper is organized as follows. In section 2 we recall the basic definitions and 
notations for multi-cut matrix models and their phase spaces. In section 3 we formulate 
the discrete string equation in terms of a generating function and derive some important 
identities. Sections 4 and 5 contain a detailed description of the large N asymptotics of 
the recurrence coefficient r n ^ in the regular and critical one-cut cases respectively. The 
classification of critical one-cut models and the connection to the Painleve I hierarchy 
contained in section 5 is relatively simple but it sets the pattern for the rather more 
complicated two-cut case in sections 6, 7 and 8. The generality of the classification 
of critical two-cut models and their connection to the Painleve II hierarchy achieved in 
section 8 requires a somewhat complicated calculation, but the final result can be stated 
concisely and is illustrated in the case of a merging of two cuts in a quartic model. The 
paper ends with a brief summary and outlook section. 



2. Multicut models 



Let us consider a model §\§ with coupling constants g. It is well known [3j HJ El |25] that 
in the limit N — > oo the support J of the density of eigenvalues p(x) is the union of a 
finite number s of real intervals (cuts) 

s 

J =\J( a j>Pj)> "i < Pi < ■ ■ ■ < a s < P„ (9) 
i=i 

where 1 < s < p = (degV)/2. The conditions that determine the actual value of s 
among those allowed by this bound can be stated in terms of a polynomial h(z) defined 
in the following way: let w(z) be the Riemann surface 



™(*)=x n(z- ( io ) 

and let w\(z) be the branch of w(z) with asymptotic behavior w\(z) ~ z s as z — > oo; 
h(z) is the polynomial part of the large z expansion of V'(z)/wi(z), i.e. 

Yl^\ = h(z) + 0(z- 1 ), z^oo. (11) 

The 2s endpoints a\, . . . , ft s in the s-cut case are solutions of the system of 2s equations 

Otj+l 

h(x)wi )+ (x)dx = 0, j — l,...,s— 1, (12) 

i 

■ VJz) 

z j^l dz = 0i J=0 ,..., S -1, (13) 

wi{z) 

h(z)w 1 (z)dz = — 47ri, (14) 



7 



where 7 is a large positively oriented loop around J. However, these equations may not 
be sufficient to determine uniquely s because they may have admissible solutions for 
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several values of s. In this case the additional condition p(x) > for x G J and the 
inequalities 

rcti 

/ h(x')w 1 (x')dx' < 0, for x < ax, (15) 

J X 

h(x')wx(x')dx' > 0, for (3j < x < <X/+i, j = 1, . . . , s — 1, (16) 

/i(i>i(i')di' > 0, for x > p„, (17) 

characterize uniquely the solution of the problem. Finally, the polynomial h(z) is related 
to the eigenvalue density by 

P(x) = ^-wx,+(x), x e J, (18) 

where Wx t +(x) denotes the boundary value of W\(z) on J from above. A model is said 
to be a regular if h(x) / on J and the inequalities ([T5l) - (fT7l) are strict. Otherwise it 
is called critical. 

As we said in the introduction we restrict our considerations to even potentials 
p 

V(X) = Y,92^, X = z 2 . (19) 

3=1 

For these models the eigenvalue support J is symmetric with respect to the origin. In 
the one-cut case J = (—a, a) and equations (fl2]) - (TH]) reduce to the single condition 



1, (20) 



A - a 2 

where we have denoted V\(X) = dy(A)/dA. Likewise, in the two-cut case J = 
(— P, —a) U (a, (3) and equations (TT2l - (fHI) reduce to 

dA V X (X) =n 

27ri v /(A-« 2 )(A-/3 2 ) ' 

dA A\/ A (A) =i 1 j 

2tt1 v /(A-« 2 )(A-/3 2 ) 

Given a family of potentials ( IT91 with the coupling constants g = (g2,gA, ■ ■ ■ ,92 P ) 
running on a certain region G of R p , we decompose the region G as 

G=\jG s , (22) 

S>1 

where g G G s if and only if g determines a s-cut regular model. The set G s is called 
the s-cut phase of the family of Hermitian models and the decomposition ( 122]) is called 
the phase diagram of the family of Hermitian models. 

To illustrate these ideas we show in figure [1] the phase diagram of the quartic 
model [51 [201 [22] 

V(X) = g 2 X + g 4 X 2 (23) 
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Figure 1. Phase diagram of the quartic model V^(A) = ^A+^A 2 in the region 54 > 0. 



in the region 

G={g=G?2,24)GR 2 :24>0}. 
The one-cut phase G\ can be written as 
G x = G^UG?, 

where 

G? = {(92,9a) eR 2 : 9 2 > 0,g 4 >0}, 
G? = {(92,94) e R 2 : g 2 < 0,24 > 0,£ 2 > -2V2l}, 

while the two-cut phase G 2 is given by 

G 2 = {(g2,9i) € R 2 : 24 > 0,22 < -2v^I}. 

The phase diagram features the critical curve 

92 = -2y/gi, 



(24) 

(25) 

(26) 
(27) 

(28) 

(29) 

which is the boundary between the two phases. 

Figure [1] also shows the oriented paths g/T traced in the phase diagram when the 
coupling constants are scaled from a given initial point g by the temperature T > 
as T — > 00. From the figure it is clear that if g e G x then g/T e Gi for all T > 0. 
However, if g e U G 2 then g/T crosses the critical curve at 



A 

424 



3. The discrete string equation 

Given N > let us denote by v n the orthogonal polynomials 

v n (x) = P n , N [x), n>0 



(30) 



(31) 
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and let us denote by L the Lax operator of the underlying Toda hierarchy 

Lv n = v n+1 + r njAr f n _i, r 0i7V = 0. (32) 

The main tool to study the asymptotics of the recurrence coefficient r n ^ is the 
discrete string equation [5] 

n 

Vz{L) n ,n-\ = Jj, (33) 

where 

p 

V z (L) = J2^g 2k L 2k - 1 (34) 
fc=i 

and the subindices (n, n — 1) denote the matrix element between the corresponding 
elements of the basis v„. We note that (13"3l can be rewritten as 



^:V x (\)U n>N (\) = ^ (35) 

27T1 JM 

where 7 is a large positively oriented circle |A| = R, and U^nW is the generating 
function 

U njN (X) = 1 + 2 ^(L 2fe_1 ) n!n _ 1 A _fe . (36) 

k>l 

In the next sections we calculate the large N expansions of the recurrence coefficient 
r n ,N by solving appropriate continuum limits of the string equation in the form fl35l) 
with a two-step procedure: (i) we first calculate a large iV expansion for the generating 
function U U) n whose coefficients are functions of the expansion coefficients for r n A r, 
and (ii) we substitute this expansion for U n> N into the continuum limit of the string 
equation (|35l) . perform the contour integration (i.e. pick up the coefficient of A" 1 in the 
large A expansion of the integrand), and solve recursively for the coefficients of the r n ^ 
expansion. 

The first step relies on two identities for U rit N(X) which we derive in the following 
theorem. 

Theorem 1 The generating function U n ^{\) satisfies the linear equation 

A(t/ ?1+ i > Af — Uu^n) = T n +l,N(Un+2,N + ^n+l,iv) — r n,N{U n ,N + Un-l,N) (37) 

and the quadratic equation 

r n ,N{U n ,N + U n _ hN )(U n:N + U n+ i >N ) = X(U^ N - 1). (38) 

Proof. For clarity in this proof we will drop the subindex N in U n ^ and denote 
the generating function ( )36l) by U n . Let us prove first the linear equation ( 1371) . In 
terms of L the recurrence relation (j6]) can be written as Lv n (x) = xv n (x). Therefore 
L 2k ~ x vAx) = x^^Vnix) and 



1 



DC 



(L fe_ ) n>n _i = / X^Vn^Vn^x)^ (39) 

n n -l,N J-oo 
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where d/i = e~ NV ( x >dx. Hence, 

u. = i + ^ f" £!^EK_iM d ^ (40) 

">n-l,N J-oo A — X 

and 

A(C/„ +1 - U n ) = - / xv n+1 (x)v n (x)d/2 - - / xv n (x)v n _x(x)dn 

iln,N J-oo ftn-l,JV J-oo 

2 Z" 00 x 3 w n+ i(x)w n (a;) 2 Z" 00 x 3 w n (x)w„_i(x) 



+ ^_ / -7^ 2 d/a - — — / -^y v M /x. (41) 

<ViV J-oo A — X IT>n-l,N J-oo A — X 

Using x j v n (x) = LPv n {x) for j — 1, 2 we find that 

xw n+ i(a;)t; n (a;)d/j = /i n+ i,jv, (42) 



oo 

00 s 3 w n+ i(x)?; n (a;) 

2 d ^ 



— oo 



A — X 

xv n+2 v n+ i + x{r n+ltN + r n , N )v n v n+ i + xr n ^ N r n _ ljN v n+ iV n ^ 2 



A -x 2 



d/i, (43) 



and 

x 3 w n (x)w n „i(x) 



A -x 2 

oo 



d/i 



L d "- (44) 

Substituting (j4"2"]) - (jl4"]) into (T4"Tj) and taking into account (T4"Uj) we conclude that the linear 
equation fl37j) holds. It is now easy to prove the quadratic equation f l38]) . Indeed, the 
linear equation f[3T|) implies 

r n +i,jv(^n+2 + f/n+i)(f/n+i + U n ) - r njN (U n+1 + U n ){U n + [/„_!) = A(£^ +1 - (45) 
and therefore the expression 

\U 2 n - r n , N (U n+1 + U n )(U n + C/„_0 (46) 
is independent of n. Since r A r = and [7 = 1 the quadratic equation (1581 follows. 

4. Regular one-cut models 

In this section we show how to calculate large N asymptotic expansions in regular 
one-cut models by passing to the continuum limit in the quadratic identity (138|) . 

4-1. Large N expansions for regular one-cut models 
Consider the continuum limit 

n oo, N ->• oo, — -> T, (47) 
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in the case where g/T G G\. It has been rigorously proved [3J [321 EZ] that under these 
assumptions there exists an asymptotic power series for the recurrence coefficient r n ^ 
of the form 

r n , N ~r(e,T) = J2n(T)e 2k , e = 1 (48) 

fc>0 

which is uniform with respect to T in a neighborhood of T = 1. Moreover, 
a 2 

r = - (49) 

where (—a, a) is the eigenvalue support for the model with coupling constants g/T. We 
write a similar series for generating function U n ,N 

U n , N (X) ~ C/(A, e;r) = J2 U k (X] r Q , . . . , r fc )e 2fc , (50) 

fc>0 

in which our notation for the coefficients U k anticipates their dependence on the rj up 
to j = k. Substituting (JlSj) . (j50j) and the corresponding shifted expansions 

r n+ii7V ~ ryj(e,T) = r(e,T + je), j G Z, (51) 
U n+hN (\)~U y] (\,e;r) = U(\,e;r b] ), j G Z, (52) 

into (135]) . the continuum limit of the quadratic identity is 

r(U + [/[_!])([/ + £/[!]) = A(f/ 2 - 1). (53) 
Identifying powers of e in (1531) we find that the coefficients Uk can be written as 



■* " V A- (84) 

= g fcy - (55) 

where Ukj is a polynomial of degree j in r , . . . , and their T derivatives. For example, 

E/i,i = 2n, C/i,2 = 2r C £/ 1>3 = 10r (r[,) 2 . (56) 

It is also easy to see that in general 

£4,i = 2r fc (57) 

and that the dependence of Uk in comes solely from Uk,i- Therefore 

U k = U (t^t- + ■■■) (58) 



\A-4r / 

where the dots stand for terms in r , . . . , rk-i and their T derivatives. 
Analogously, the continuum limit of the string equation (T35|) reads 



^V x (\)U(\,e,r) = T (59) 

Z7T1 
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or using the expansion (I50|) for U, 

{ |Va(A)£4(A; r , . . . , r k ) = 8 kfi T, k > 0. (60) 

J 7 27T1 

The first of these equations (corresponding to = 0) can be written as 

W(r Q ) = T, (61) 

where 

w ^-i^^lSrtX") kS2kr °- (62) 

Equation (ISTI) is an implicit algebraic equation for r of hodograph type (linear in the 
variables g and T) which is equivalent to the normalization condition ( 120]) due to the 
identification ()49l) . 

The remaining equations fl60|) can be written as 

3k 

^ Cj{r Q )U k j{ro, . . . , r fc ) = 0, fc > 1, (63) 
i=i 

where 

and where W^(r ) = d j W(r )/dr J . 

For an even potential in the one-cut case h(z) is an even function of z whose 
characterization ffTTj) in terms of A = z 2 reads 

n VaVa(A) 



2 



/i(A) + 0(A' 1 ), A^oo, (65) 



and consequently 



Wto) = tf£ .f™ =h(X) . (66) 

/ 7 2vri (A - 4r ) 3 / 2 A= a 2=4r 

Thus, if g/T G Gi then W'(r Q ) ^ and the implicit function theorem shows that flBTT) 
determines locally smooth function of T. Moreover, from the dependency of U k 

on rfc shown in (lo^j) it follows that equations (1631) are of the form r k W'(r ) = • • •, where 
the dots stand for a sum of terms on r , . . . , 7"fc_i and their T derivatives. Therefore 
the implicit function theorem can also be applied and ( 163]) determine recursively all 
the coefficients r k of the large N expansion ( HHl) as locally smooth functions of T. 
Furthermore, if we differentiate the hodograph equation (16T|) with respect to T we can 
write the T derivatives of r k as rational functions of the W^' and hence as rational 
functions of r and g. For example, from r' = 1/W'(r ) and the k = 1 equation (163|) 
we find 

= / W'\r f _ ^ w (r ) \ 
ri V6^'(r ) 4 12^'(r )V r °' 1 J 
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4-2. The quartic potential in the regular one-cut case 

As our first application of these results consider the quartic potential ( 1231) in the regular 
one-cut region. Equation ( 1621) gives 

W(r ) = 2g 2 r + 12(? 4 r 2 . (68) 

The leading term ro(T) of the asymptotic expansion (1481) is the positive root of the 
hodograph equation W(r ) = T, namely 



r„(T) = (69) 

To calculate the next term r\ (T) we first use (!67|) and (1681) to obtain r\ as a function of 
7o, and then we substitute the former explicit expression of r$(T) to get 

( rj* g^g 2 2 + l2Tg A - g 2 ) 

ri(T) = 2 {3 i + i2T 9l y- ■ (70 » 

Higher coefficients r k {T) of the asymptotic expansion ( 148]) can be easily calculated 
iterating this procedure. We postpone a graphical illustration of r (T) until we discuss 
the two-cut region of the same model in section 6. 

4-3. The Brezin-Marinari-Parisi potential in the regular one-cut case 

As our second example we consider the sixtic potential 

V(X) = 90A - 15A 2 + A 3 , (71) 

introduced in [28J by Brezin, Marinari and Parisi to generate a non-perturbative 
ambiguity-free solution of a string model. The matrix model corresponding to V(X)/T 
is regular except at T c = 60, but this critical temperature does not mark the boundary 
between a one-cut region and a two-cut region. The function W takes the form 

W(r ) = 180r - 180r 2 + 60r 3 , (72) 

and the unique positive root of the hodograph equation W(r ) = T is given by 

r (T) = 1 + S/T/60 - 1. (73) 



r l( T ) = eJnXJen TvP ( 74 ) 



Proceeding as in the previous example, equations (1571) and ( 1731) yield 

1 + {/T/60 - 1 
64800(T/60 - 1) 

In figure [2] we plot r (T) as a function of T in a neighborhood of the critical T c = 60. 
Note how the vertical tangent to the graph marks a critical point which is qualitatively 
different from the graph of the quartic potential that we will discuss in section 6. 
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Figure 2. Leading coefficient ro(T) as a function of T for the Brezin-Marinari-Parisi 
potential V(A) = 90A - 15A 2 + A 3 with a critical point at T c — 60. 



5. Critical one-cut models and the Painleve I hierarchy 

In this section we first give a precise definition of singular solution of order m in a one- 
cut model, then we show how the algorithm to calculate the asymptotic expansion of 
the recurrence coefficient developed in the previous section has to be modified by means 
of a suitable double scaling, and finally we give a general proof of the relation between 
a one-cut singular model of order m and the (m — l)-th Painleve I equation. 

5.1. Critical one-cut models and double scaling limit 

Definition 1 A solution (j-q, T) = (r c , T c ) of the hodograph equation W(ro) = T is 
called a singular solution of order m > 2 if 

W(r c ) = T c , W\r c ) = ■■■ = W {m ^ (r c ) = 0, W ^ (r c ) ^ 0. (75) 
The identity 

H*> W = (2k- 1)112* i± ^ , = ^I^T'^tMA) (76) 
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shows that singular solutions of order m of the hodograph equation correspond to critical 
matrix models for which the function h(z) has zeros of order m — 1 at the endpoints 
of the eigenvalue support J = (—a, a). We note also that the eigenvalue density p(x) 
verifies 

p[x) ~ (x=Fa) m "5, z-)>±a. (77) 

These critical matrix models correspond to conformal (2m— 1, 2) minimal models [7]IT6]. 

Let us thus consider a model that (locally) for T < T c is in the regular one-cut case 
and that for the critical temperature T = T c the corresponding hodograph equation (I6T1) 
has a singular solution r c of order m > 2. Since W(r c ) = we cannot invoke the 
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implicit function theorem to solve ( 16 ip near T = T c with r (T c ) = r c . In fact, r (T) can 
be expanded in powers of (T — T c ) l l m and therefore the system (160]) does not determine 
the coefficients of the expansion ( 1481) as locally smooth functions of T at T = T c . We 
regularize this critical behavior introducing a scaling variable x and a scaled expansion 
parameter e defined by 

- _ e l/(2m+l) _ jy-l/(2m+l) 5 (73) 

T = T C + e 2m x. (79) 
In terms of these scaled variables the string equation ( 1591 reads 

/ ^V X (\)U(\, e; r) = T c + e 2m x, (80) 

/ 7 2l7T 

and we look for an asymptotic expansion of the form 

r(e,x) = r c + ^x k {x)t 2k . (81) 
fc>i 

The shifts T — >• T ± e correspond to x — > x ± e and therefore U(X, e; r) is determined by 
the quadratic equation 

r(U + U FT] )(U + U [I] ) = \(U 2 -l), (82) 

where we are denoting fm(x) = f(x + ke). 

The rest of the calculation is straightforward and can be carried out in complete 
analogy with the previous section. The corresponding expansion of U is 

U(X,e;x) = J2u [k] (Kx 1 ,...,t k )e 2k . (83) 

fc>0 

Substituting ( IHTi) and ( 1831) in ( 1821) and equating powers of e we find that 



y w = ^r yM ^---'^i) fc>1 (85) 

{\-Ar c y - 

and that the coefficients U^ k ^ are polynomials in ti, . . . , tk-j+i an d their x derivatives, 
which can be determined recursively [23] • In particular 

U [k ' 1] = 2x k . (86) 

The first few arc 

U [2 ' 2] = Qxl + 2r B x?l, (87) 

jj[3,2] = 12rit2 + 2r ir " + 2r c x" + -r c tj 4) , (88) 

6 

^[3,3] = 2Qr 3 + lOr^r;) 2 + 20r c titi + 2r 2 x { ^ . (89) 

In the following theorem we prove that the diagonal coefficients U^ k ' k ^(xi) are the well- 
known Gel'fand-Dikii differential polynomials of the KdV theory and that the coefficient 
Xi(x) is a solution of a member of the Painleve I equation. 
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Theorem 2 The coefficient u = X\ for a singular one-cut model of order m satisfies the 
(m — \)-th Painleve I equation (Pj l ~ l equation) 

c m (r c )R m (u) = x, (90) 

where R m (u) is the differential polynomial in u determined recursively by 

d x R m+1 = {r c dl + Aud x + 2u x )R m , R = 1, (91) 

and 

Cm (r c ) = -11 (92) 

v ' 2 m (2m- 1)!! v ; 

Proof. The linear equation 

Tft{U m + U m ) - r(U + Up^) = X(U {1] - U) (93) 

follows immediately from the quadratic equation (152]). and substituting the expansion 
for U given by (jS3]l and f[55"j) into (JS3J) we find that the diagonal coefficients U^- k,k ^ satisfy 
the recursion relation of the Gel'fand-Dikii differential polynomials, 

d jjlk+iM-i] = ( rc5 3 + AxiQx + 2 ( Xl ) x )U [k ' k \ U m = 1. (94) 



If we now substitute the expansion (l83l) for Z7 into the scaled string equation (1801 . equate 
powers of e, and take into account f l64"|) as well as the definition of singular solution f J75|) . 
we obtain 

fc 

c^rc)^^ 1 ^!, . . . , x k - j+ i) = 5 kym x, k>m, (95) 



so that 



c m (r c )f/[ m ' m ](n)=a;, (96) 



m+k 



CjtrjU^Wfa, . . . , r m+fc _ i+1 ) = 0, k > 1. (97) 



3=m 



This system provides for each coefficient tfc(x) (A; > 1) an ordinary differential equation 
involving the previous coefficients tj-(l < j < k). In particular (|96|) takes the form (1901) 
as stated. The set of ordinary differential equations (1901 constitutes the Painleve I 
hierarchy, whose m = 2 member is the familiar Painleve I equation 

c 2 {r c )(2r c u xx + 6u 2 ) = x. (98) 

Note that a more precise notation for R m (u) and c m (r c ) in theorem 2 would be 
R m (r c ,u) and c m (r c ,g), which emphasizes the explicit dependence on the critical value 
r c and the coupling constants g. 
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5.2. The critical Brezin-Marinari-Parisi model 

We can apply immediately Theorem 2 to the Brezin-Marinari-Parisi potential. The 
singular solution (r c = 1,T C = 60) of the corresponding hodograph equation is of third 
order, since 

W(r c ) = T c , W'(r c ) = 0, W"(r c ) = 0, W"\r c ) = 360. (99) 
Therefore u = ti is as solution of the Pj equation 

u xxxx + I0uu xx + hu 2 + 10m 3 = -x. (100) 

6 

6. Regular two-cut models 

In this section we study regular two-cut models (i.e. g/T G G 2 for a given T > 0) 
following the same pattern used in our previous study of regular one-cut models. The 
main difference is that now, due to ([7]), we need two distinct asymptotic expansions: one 
expansion for the odd recurrence coefficients r 2n +i,N and a second expansion for the even 
recurrence coefficients r 2n ,N- This fact leads us to introduce two different generating 
functions V and W for U n ^ with n odd and even, respectively. It also requires the 
splitting of the quadratic equation ( )38l) for U n at into a system of two coupled equations. 



6.1. Large N expansions for regular two-cut models 

In view of ([7]) we formulate the asymptotics of r n ^ in the limit 

2n 

n^oo, N 00, Jj^ T ( 101 ) 

as 

r<2n+i,N ~ a[i](e,T), r 2n ,N ~ b(e,T), (102) 

where we use again the bracket notation for shifts and we assume that a(e, T) and b(e, T) 
are asymptotic power series 

a(e,T) = ^a fc (T) e 2fe , b(e, T) = £ b k (T)e 2k , e =l (103) 

fc>0 k>0 

We notice that contrarily to the one-cut case and with the exception of the quartic 
model [5j [20J, [22], to our knowledge there is no rigorous proof that power series with 
smooth coefficients are truly asymptotic to the recurrence coefficients. 

To calculate the coefficients a k and b k we perform the continuum limit fllOip in the 
string and quadratic equations ( 135]) and ( |38|) . and introduce two generating functions 

2k 



V(A, e; a,b) = V fc(^ a o> b o, ■ ■ ■ , «fc, h)^ 



k> ~" (104) 
W(A, e; a,b) = w fcO; a o, K ■ ■ ■ , b k )e 2k , 

k>0 
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such that 

£W,iv(A) ~ V [1] (A,e;a,6), U 2n , N (\) ~ W(A, e; a, b), (105) 
and more generally 

t/ 2 (n+i)+i,7v(A) ~ V[ 2 j+i](A,e;a,6), t/ 2 (n+j),iv(A) ~ W [2i ](A, e; a, b). (106) 
Then (|38|) is equivalent to the system 

«(V + W M] )(V + W W ) = A(V 2 - 1), 

6(W + V h i])(W + V[i]) = A(W 2 - 1). 
Identifying the coefficients of e° in these equations we get 

a - b + A 6 - a + A 

V = , W = , (108) 

w w 

where 



w = VA 2 - 2A(a + 6o) + (&o - ao) 2 = \/(A - « 2 )(A - /3 2 ). (109) 

The last equality follows from ([7]) , which allows us to express Oq and &o i n terms of the 
endpoints of the eigenvalue support J = (— (3, —a) U (a, (3) of the model with coupling 
constants g/T: 

a, = - A (a + (3)\ 6 = l(a-/3) 2 . (110) 

Recursive identification of the coefficient of e k for k > in f ll07p leads to a system 
of two linear equations in and of the form 

(AV - 2w~ 1 \a )V k - 2w- 1 Xa W k = ■■■ 

-2w- l \b Q \ k + (AW - 2w- 1 \b )W k = ■■■ 

where the right-hand sides are functions of ao, b , . . . , a k , b k and their T derivatives. 
Since the determinant of the coefficients of (11 111) equals A 2 we can solve uniquely (II 111) 
for V k and as functions of a ,b , . . . ,a k ,b k and their T derivatives. Moreover, if 
V(A, e; a, b) and W(A, e; a, b) satisfy (llOTj) . so do the functions 

V(A,e;a,6) = W(A,e;6,o), W(A, e, a, b) = V(A, e; b, a). (112) 
From these equations and taking into account that as A — > oo 

V(A,e;a,6) ~ 1 + 0(A" 1 ), W(A, e; a, b) ~ 1 + 0(A" 1 ), (113) 
it follows that V and W are related by 

W(A,e;o,6) =V(A,e;6,o). (114) 

Using induction in the quadratic system (11071) we find that the coefficients can 
be written in the form 

v* = £^t^< (115) 

3=0 
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where R kJ = R k ,j(a , b , . . . , a k , b k ) and S kJ = S ktj (a , b ,..., a k , b k ) are polynomials in 
ao, bo, ... , a k , b k and their T derivatives. It is immediate to prove that the dependence 
of V fc on (a k , b k ) is given by 



2a k (2A(a + bo) - 2(b - a ) 2 )a k + 4Aa &. 



V k = —+ ^ u ' - ' "-""' fc + ■ ■ ■ , k>\, (116) 



w w 3 



where the dots stand for terms dependent on a , b , . . . , a k ^%, b k _i and their T derivatives. 

The next step is to calculate the a k and b k as functions of T. Substituting the 
expressions for V and W into the continuum limit of the string equation ( 135]) it splits 
into the system of equations 



j^V x {X)V{\,e;a,b) = T, 
/^(A)V(A,e;M = T, 



(117) 



or, equivalently, 



7r4vx(A)V fc (A; oq, b ,..., a k , b k ) = 5 kfi T, 

Z7T1 

^:Vx(A)V fc (A; b , a , ■ ■ ■ , b k , a k ) = 5 kfi T. 

Z7T1 



(118) 



For k = we have a system of two hodograph equations for do(T) and bo(T), 

W(a ,b ) = T, W(b ,a )=T, (119) 

where 

W(a ,b )= /^(A) A + Qo ~ & ° . (120) 

J 7 27T1 W 

Note that in view of (IllOp the equations (I119p are equivalent to the equations (j2ip which 
determine the two-cut support for a model with coupling constants g/T. 

Finally, we say that (do, bo) is a regular solution of f lll9p if the Jacobian determinant 
of (I119P with respect to (cto, b ) does not vanish. Otherwise (do, bo) is called a singular 
solution. Equation (I116p shows that if (d , b ) is a regular solution of f )119p . then the 
system f II 1 8 j) determines recursively all the coefficients a k and b k as locally smooth 
functions of T. 



6.2. The quartic model in the regular two-cut case 

As we have discussed at the end of section 2, for g 2 < and T < T c = g^/i^g*) the 
quartic model (1231) is in the regular two-cut case, and we can calculate expansions of 
the form fl 1 3 [) . The first coefficient Vo can be read off directly from (I108p . Thus, the 
string equations (11181) for k = are 

2# 2 ao + 4^ 4 (do + 26 d ) = T, 
2g 2 bo + Ag 4 (bl + 2a bo)=T, 
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and we find immediately the leading coefficients 



ao( r) = Va-ffr-a (122) 

bo (T) = -yfr 2 -*^-^ 

4#4 

Then we solve f llOTj) for the next term Vi, of which in the notation of equation (11151) 
we need the following coefficients: 

#1,0 = 2ai, 5i i0 = 0, Si A = 2(a + b )a 1 + 4a &i + 2a b' \ (124) 
With these expressions we find that the string equations (11181) for k = 1 are 

2g 2 ax + 2g 4 (4a ai + Ab ai + 4a &i + 2a & ') = 0, 

2# 2 &i + 2^ 4 (46 &i + 4a 6i + 46 ai + 2b a'^) = 0. 
Hence, the second coefficients are 



(rr* gM + 4Tg±-g 2 ^/g 2 2 -4:Tg 4 ) 
fll(T) = 2{g 2 -ATg 4 fl 2 > ( 126 ) 

HT) = 2(rf-4W • (127) 

In figure [3] we illustrate these asymptotic behaviors for a quartic model with g 2 = —2 and 
#4 = 1, so that the critical temperature is T c = 1. For < T < T c we plot the leading 
coefficients ao(T) and bo(T) of the asymptotic expansions for the odd and even terms 
of the recursion coefficient in the two-cut region given by equations d 1 2 2 [) and (11231) 
respectively, while for T > T c we plot the leading coefficient tq(T) of the asymptotic 
expansion in the one-cut region given by equation ( 1691) . As the temperature increases 
through T c the model features a merging of two cuts [20]. In fact this is the simplest 
example of the general processes that we classify in the next section. 

7. Critical two-cut models with merging of two cuts 

In order to classify critical models in the two-cut case it is useful to consider the function 

f d A 

F{ao, b ) = & —V x (X)w(X, ao, b Q ) + T(a + b ), (128) 

J 7 27T1 

where 

w(X, oq, b ) 2 = A 2 - 2A(a + 6o) + (bo - a ) 2 , (129) 
because the solutions of the system 

dF _ o dF - o (130) 
da ' <% 

are precisely the solutions of the hodograph system (11191) . Moreover, it is convenient to 
change variables from (do, fro) to a new pair (a, t) defined by 



o = a 2 = a + b - 2^/a b , r = (5 2 = a Q + b + 2^/a bo, (131) 
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Figure 3. Leading terms ro, ao and 60 as a function of T for the quartic model 
V(A) = .92^ + ff4^ 2 with 52 = — 2 and 34 = 1. 



where a and (3 determine the endpoints of the eigenvalue support J (cf. (IllOp ). In these 
new variables the function ( 11281) is given by 

F(a, r) = j> ^ [ V x (X) v /(X-a)(X-r) + T -{<j + r), (132) 

and satisfies the Euler-Poisson-Darboux equation [29] 

d 2 F dF dF 

= < 133) 

Hence it is clearly advantageous to analyze the solutions (ao,&o) of the hodograph 
system (11301) in terms of the solutions (a, r) of the transformed system 

dF dF . , 

a7 = °- *=°- < 134 » 

because for any solution (a, r) of (11341) with a 7^ r the Euler-Poisson-Darboux equation 
permits to express any mixed derivative d l+ ^ F jda % dr^ as a linear combination of 
unmixed derivatives d n F/da n and d m Fjdr m . 

7.1. Merging of two cuts 

We are interested in solutions of (11341) with a = (i.e. ao = bo). Since 
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we have that 



d 2 F 
dal 



d 2 F 



d 2 F 
da db 

d 2 F 



db 2 



d 2 F , OF . 
32^(0, r)^(0,r) 



0. 



(137) 



da db 

Therefore these solutions of (I134p with o = determine singular solutions of (I130p . 
Moreover, the corresponding matrix models have an eigenvalue support with a = so 
that they represent a critical process of merging of two cuts. This motivates our next 
definition: 



Definition 2 A solution (a, r) of Ijl34\ ) with a = and r > determines a singular 
solution of the system ( fiffflj) with a merging of two cuts of order m > 1 if 

Qk p Qm+l p 

.(0,r) = 0, fc = l,...,m; ^-^(0,-^0, 
d 2 F 



dF 



(138) 



(0,r) = 0, 



-(0,r)^0. 



dr ' ' Or 2 

For these solutions of ( 11301) we have that 

(2k- 1)!! r dA (A — r)V\ 



ak+l p 



2vri A fc 



w 



(2Jfe-l)!! d 



2 fc + 1 (fc - 1)! SA^- 1 
where /i is the polynomial in A = z 2 given by 



(A - r)h(X) 



A=0 



h(A) 



2z 



k > 1, (139) 



(140) 



Thus /i(^) has a zero of order 2m — 1 at z = and the eigenvalue density p(x) verifies 

p(x) ~ x 2m , a; 0. (141) 

These critical models are related to conformal (2m, 1) minimal models |16[ IT?] . 

Note that there is a difference between the definitions of order of a singular solution 
in the one-cut (definition 1) and two-cut (definition 2) cases, which in turn entails a 
difference in the statements of our main results (theorems 2 and 3). However, with 
these definitions we achieve a complete analogy in the form of both scaling expansion 
parameters e (ITS]) and f l 1 5 1 j) as well as the corresponding string equations (180]) and (j!60p . 



7.2. Merging of two cuts in the quartic model 

As an illustration of these ideas, consider the quartic model (|23|) . The function ( I132p is 

F(a, r) = ^(-g^o- 3 - g 2 a 2 + 9i ra 2 + 2g 2 ra + g 4 r 2 a - g 2 r 2 - 9i T 3 ) + ~T(a + r). (142) 

It follows at once that at the critical temperature T c = g\j{^g±) we have a solution 
(0, r c ) of the system 



OF 
'do 7 



(0. 



Tr 



0. 



OF 



(0,T C 



0. 



(143) 
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given by r c = — g^/ffe- Moreover, since 

|^(0, t c ) = -f ^ 0, (144) 
this solution describes a merging of two cuts of order m = 1 . 

7.5. Double-scaling limit and string equations 

Let us consider a matrix model V(X)/T with a critical point at T = T c such that locally 
for T < T c the model is a regular two-cut model, and at T = T c the model features a 
singular two-cut merging of order m as described in the previous subsection. Then 

T^-^^T^'J^^r- < 145 » 

so that the two cuts of the eigenvalue support merge to J = (—f3 c , (3 C ) where (3% = 4r c , 
and the hodograph system (I119p reduces to 

|Va(A)-=T c , 

w c 

dX 1 (146) 

Z7T1 W r 



where 



w c = v / A(A-4r c ). (147) 

Note that the hodograph system (11461) in effect determines both T c and r c in terms of 
the coupling parameters g. Moreover, the definition (11381) of singular solution with a 
merging of order m reads 

dA (A - 4r c ) k+1 V x 



7 27ri wl k+l 



0, k = l,...,m- 1, (148) 



7 27ri w c 2m+1 ' / 7 2vri w% 

Equations (I148j) and (11491) characterize the subset of the phase space G representing 
these critical models. 

As in our study of critical one-cut models, the implicit function theorem cannot be 
applied to solve ( 159"|) near T = T c , with a (T c ) = b (T c ) = r c . In analogy with (179"]) we 
regularize this critical behavior in the continuum limit 

2n 

n^oo, iV->oo, J^^ T ( 15 °) 
using a scaling variable x and a scaled expansion parameter e defined by 

- _ e l/(2m+l) _ N -l/(2m+l)^ 

T = T C + e 2m x. (152) 
The continuum limit of the recurrence coefficient can be written as 

r 2 n+i,N ~ 0|i|(e,x), r 2 „ >A r ~ b(e,x), (153) 
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where we denote fm(x) = f(x + ke) and a(e, x) and b(e, x) denote asymptotic power 
series 

a(e,x) =r c + ^2a k (x)e k , b(e, x) = r c + ^ b k {x)e k . (154) 

k>l k>l 

We introduce again two generating functions 

V(A, e; a,b) = J2 a 1; b 1; . . . , a fc , b k )e k , 

^ (155) 
W(A, e; a, b) = ^ W W (A; a 1} b 1; . . . , a fc , b fc )e fc , 

fc>0 

such that the asymptotic behaviors of the odd and even terms of the generating function 
U n are given respectively by 

EWi(A)~ V m (A,e;a,b), f/ 2n (A) ~ W(A, e; a, b). (156) 
Thus the continuum limit of the quadratic equation ([38]) for [/„ splits into the system 

a(V + WpqJCV + W m ) = A(V 2 - 1), 

b(W + Vp^)(W + V m ) = A(W 2 - 1). 
Identification of the coefficients of e° in fl 1 5 T[) leads to 

V [0] = W [0] = A, (158) 

and recursive identification of the coefficients of e k for k > to a system of two linear 
equations for and of the form 

(2A - 4r c )V^ - 4r c W^ = ■ • • , 

_4r c V [fc] + (2 A - 4r c )W [fc] = • • • , 1 °' ' 

where the right-hand sides are functions of Oi, bi, . . • , Ofe, b^ and their x derivatives. Once 
more, the determinant of the coefficients does not vanish and (11571) uniquely determine 
the coefficients and W' fc ' as polynomials in the functions ai, bi, . . . , a&, b k and their 
x derivatives. 

Finally, the continuum limit of the string equation ( l35l) splits into the system 



j> ^^a(A)V(A, e; a, b) = T c + e 2m x, 
/ ^^(A)W(A, e; a, b) = T c + e 2m x. 



(160) 



8. Critical models with merging of two cuts and the Painleve II hierarchy 

In this section we finally show how a particular class of solutions of the string equations 
for two-cut merging models is related to the Painleve II hierarchy, but before embarking 
on this calculation we anticipate informally our result, recall briefly our two-step 
procedure to calculate the formal asymptotic expansion of the recurrence coefficient 
and discuss the source of the technical complications we will have to deal with. 
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Theorem 3 will essentially state that the coefficient ai(x) in the formal asymptotic 
expansion of the recurrence coefficient (I154p for a symmetric solution of a two-cut 
merging model of order m satisfies the m-th Painleve II equation. The general 
strategy to prove this theorem is again (i) to calculate large N expansions for the 
generating functions V and W whose coefficients are functions of the and b k (and 
their derivatives) by using the continuum limit of the "separated odd-n even-n" version 
of theorem 1 given by equations (I157p . and (ii) to substitute these expansions for V 
and W into the continuum limit of the "separated odd-n even-n" version of the string 
equation (1351) given by equations (11601) . perform the contour integrations, and solve 
recursively to obtain differential equations for the coefficients a k and b k of the odd-n 
and even-n terms of the r n ^ expansions. This recursive solution is, however, difficult to 
carry out in full generality because the order k in ( 11551) at which the x dependent term 
in (11601) enters the expansions depends on the order of the critical model. 

8.1. Symmetric solutions of the string equations for two-cut merging models 

In this section we discuss asymptotic expansions of the recurrence coefficient (11541) for 
which b k (x) = (— l) k a k (x) (symmetric solutions), and show that for these solutions the 
odd and even generating functions V and W can be replaced by a single generating 
function V. 

Definition 3 A solution a(e,x), b(e, x) of the system $16(A) is said to be symmetric if 

b(e, x) = a(— e, x). (161) 

Note that if the functions V(A, e; o(e, x), a(— e, x)) and W(A, e; a(e, x), a(— e, x)) satisfy 
the scaled quadratic system (11571) so do the functions 

V(A, e; a(e, x), a(-e, x)) = W(A, -e; a(-e, x), a(e, x)) 

_ _ (162) 

W(A, e; a(e, x), a(-e, x), ) = V(A, -e; a(-e, x), a(e, x)). 

Hence symmetric solutions satisfy 

W(A, e; a(e, x), a(-e, x)) = V(A, -e; a(-e, x), a(e, x)). (163) 
Consequently, it is convenient to introduce the function 

V(A, e; x) = V(A, e; o(e, x), a(-e, x)) (164) 
since the scaled string equations (I160p reduce to a single equation for V, namely 

j ^V x {\)N [h \\- tlx,..., Ok) = 6 kfi T c + 5 k>2m x, k>0. (165) 

The function V has an expansion of the form 

V(A, e; x) = ^ ^ (A; a u . . . , a k )e k , (166) 

fc>0 

and the corresponding coefficients are determined using the first equation of (11571) 
a(e, x)(V(A, e; x) +V(A, -e; x- e))(V(A, e; x) + V(A, -e; x + e)) = A(V(A, e; xf - 1).(167) 



y[0] = 


\ 

A 






V [i] = 


2ai 




w c 


V [2] = 


2a 2 




w c 


V [3] = 


203 
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For example, the first few coefficients are 

(168) 
(169) 

+ A ^ a V 2 ^ (170) 

wl 

| 8rX-4ag + A(4a 1 a 2 -2r c a?) ^ 

w c wl 

8.2. Structure of the coefficients 

To proceed further we need a general expression for the structure of the coefficients 
V[ fc l as functions of Oi, . . . , Ofe and their x derivatives, which could be derived by a 
direct analysis of the quadratic equation (11671) . Although straightforward in principle, 
in practice the resulting intermediate expressions are complicated. We achieve a certain 
simplification by forming suitable linear combinations of shifted equations that exhibit 
well-defined parity. In fact, we resort to this method twice (propositions 1 and 2). 

Proposition 1 The following linear expressions in V are even functions of e: 

AV(A, e; x - e/2) - a(e, x - e/2) [V(A, e; x - e/2) + V(A, -e; x - 3e/2)] , (172) 
AV(A, e; x + e/2) - a(e, a; + e/2) [V(A, e; a; + e/2) + V(A, -e; 2 + 3e/2)] . (173) 

Proof. If we perform the shift x — > x — e/2 in ( I167p we get 

A(V(e; x - e/2) 2 - 1) = a(e, x - e/2) [V(e; x - e/2) + V(-e; x - 3e/2)] 

x [V(e; x - e/2) + V(-e; x + e/2)] . (174) 

The difference between this equation and its version with the substitution e — > — e yields 

A [V(e; x - e/2) - V(-e; +e/2)] = a(e, x - e/2) [V(e; x - e/2) + V(-e; x - 3e/2)] 

- o(-e, x + e/2) [V(-e; x + e/2) + V(e; x + 3e/2)] (175) 

which means that (I172p holds. Similarly, if we introduce the shift x — > x + e/2 in (11671) 
and perform the difference between the resulting equation and its version with the 
substitution e — > — e we get (I173p . 

The vanishing of the coefficients for odd powers of e in (11721) and (11731) provide 
us with a series of A-dependent constraints. To take advantage of these constraints we 
first make explicit the A dependence of V in a convenient form. Thus, using recursion 
in (11671) we deduce that the functions V^' can be written as 

% > 1, 




WC V 3=1 

= — (c^ + (fjW^T +1] + 9j(X)B : 
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where 



fj( X ) = n = — \i ' 9jW 



A 



w 2j M ' yJV ' w 2 J (A-4r c )i 



25 



(177) 



and Af ] ,Bf ] ,CW are A-independent polynomials in oi, . . . , Ofc and their x derivatives. 
Then V can be expressed in the form 



V=-(A + V + X;(/i(W 

\ J>1 



?+^(A)V?' ] 



where 



V = ]T ^e fe , V M = Af j+i] e 2j+i , Vf = ]T B f j+i] e 2j+i 



:i78) 



(179) 



fe>i 



i>0 



i>0 



For example, from fjlTOj) and f 1 1 7 1 j) we find 



Af 1 = 0, 



B [2] = ^-(4r c a 2 -ai), B\ 



A[ 



2 Ul ArV 



[3] 



4r;? 



(4r c a 2 - a?) 



180) 



r [2] _ 

" 2r ' 



C m = 2a 3 -±a'i + ^(a 2 -2r c a 2 ) 



Proposition 2 The following expressions are even functions oft: 
V (e; x - e/2) + V^e; x - e/2) + V? ] (e; x - e/2) - 2a(e, x - e/2), 
4r c v!, 1] (e; x - e/2) + a(e, x - e/2) [v (e; x - e/2) + V (-e; x - 3e/2) 



(e; x - e/2) - a(e, x - e/2) vjfte; x - e/2) + V^ J (-e; x - 3e/2) 



0], 



:isi^ 



4r c vf ] (e; x - e/2) + V? +1] (e; x - e/2) 

- a(e, x - e/2) M jl (e; x - e/2) + vf ] (-e; x - 3e/2) 

Proof. It is enough to substitute (11781) into the expression (11721) and identify coefficients 
in the A-dependent functions A, 1, fj(X), gj(X) (j > 1) taking into account that 

A/i = A-4r c , \g 1 = \ + 4r c g h 

X fj = fj-i, X 9j = -4r c gj + gj-i, j > 2. 

Equating to zero the coefficients of the odd powers of e in the third and fourth 
expressions (11811) we obtain a series of equations involving the functions A^\ and 
Oj(i = 1, . . . , k). Moreover, since f!173j) follows from (11721) under the substitutions 

V(e; x) -> V(-e; x), a(e, x) -> a(-e, x), (182) 

then the relations provided by (11731) are those supplied by ( 118ip with the substitutions 

a, t ->• (-1)^. (183) 



Af -> (-l)M/ 
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Setting to zero the coefficients of e 2 i +l and e 2j + 3 in ( |18ip we obtain 

- ai Af 3] +r c (Af j] y = 0, (184) 

-2a,Af +2] + A^ ] ~ *i{Af + 2r c (A l r 2] )' ~ + r c {Af +1] f = 0, (185) 

-a x Bf j] + 2r c Bf j+1] = 0, (186) 
-2o 3 5f 1 - 2 fll sf +2] + 4r c 5f +3] + + sf ] a 2 

+ 2a 2 (I?f ] )' - ai(sfl)' - ax(sf +1] )' - (187) 



2 

- ififa? - |ax(Bfl )" + |r c (Bf +1] )" + ^ c (^)'" = 0. (188) 

If we now sum and subtract these equations with their corresponding versions under the 
substitution ( I183p we get for j > 1 

Af j 1 = 0, (189) 

2r c cUf +2] = ai £Uf (190) 

4^ = -r c d 2 x Af +1] + 2Mf +5fl , (191) 

2r c Bf +1] = aifif ] , (192) 

2r c d x Bfi+ 2] = (2r 2 c d 3 x + 2(4r c a 2 - a 2 )^ + (4r c a 2 - cr),)/^'. (193) 

Furthermore, using the expression (TT8Q|) for .Bj 21 we have that fTT93|) can be rewritten as 

= (r c dl + 25^ + (M 2l ),)4 2il . (194) 

8. 3. Symmetric solutions of the string equations and the Painleve II hierarchy 

Let us see now how the Painleve II hierarchy emerges from the string equation of two-cut 
merging models. Substituting (1178ft into (1165ft and identifying powers of e we obtain an 
infinite series Ej (i > 0) of systems of two equations. The system E is ( I146p . and the 
system E< for i > 1 is given by 

j 

^(<^Ap(ai, . . . , om) + 7jBf\ai, a 2i )) = 5 ijm x, 

j T (195) 

J>^f +1] ( ai) . . . , a 2m ) + +1] (ai, • • • , a 2i+1 )) = 0, » > 1, 
i=i 

where 

«=i&< A »- J£L (196) 

Moreover, in view of the constraints ( 1148ft - (1149ft for two-cut merging singular models of 
order m we have 

cp 1 = ... = V 9 m _ 1 = 0, ip m ^ 0, 7i 7^ 0. (197) 
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Therefore the first m systems (I195p reduce to 



i> B f +ii = 

for i — 1, . . . , m — 1, and 



(198) 



m 

} [2m+l] 



(199) 



^r +ii +E^r +1 -°' 

for i = m. We are now ready to prove the following theorem: 

Theorem 3 The coefficient u = <Xi for a symmetric solution of a two-cut merging model 
of order m satisfies the m-th Painleve II equation (Pf} equation) 

2r c <p m (r c )R m (u) + xu = 0, (200) 

where -R m (w) is the differential polynomial in u determined recursively by 

2r c d x S m uO x R mj 



(201) 



with 



and 



= -t^u*, So = ~^u\ (202) 



2r r 8r 



c 



^ c) ^ 7 M^p • (2 ° 3) 

Proof From (11 98ft we see that the system Ei is 

7iM 2] =0, 7iM 3l = 0, (204) 
so that B[ = B[ = . Then using f ll 901) and (11941) recursively we get 

B [2i] = B [2 j+ i] = 0; j > x. (205) 

Equating to zero the coefficients of e 2j+5 and e 2 ^ +1 in the fourth expression of ( I18ip we 
obtain 

2r c Bf +3] = fll sf +2] , (206) 

= (-r c ^ 3 + 2B?d x + (B?) x )Bf^. (207) 
Then, since the system £2 is 

7 iM 4] = 0, 7iM 51 = 0, (208) 
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it follows that Bf ] = Bf = 0. Hence, using (12061) and ( 12071) recursively we get 

B f +2] = Bf +3] = 0, j>l. (209) 
Repeating this process using the systems up to i = m — 1 we get 

B [2m] = B [2m+1) = ^ j = 2 , . . . ,m, (210) 

and 

2r c Bf m+1] = a x Bf m] . (211) 
Moreover, taking into account (11891) . the system E m reduces to 



<PmA^ + 7i5l 2m+1] = 0. 



(212) 



Therefore, using (12111) we get 

2r c ip m A [ l m+l \*i) + aix = 0, (213) 

which in view of (11901) and f 1 1 9 1 j) proves that a\ is a solution of the (m — l)-th member 
of the Painleve II hierarchy ( 12011) with 

R m (u)=At m+1] (u), S m = Al m+2 \u). (214) 

Moreover, taking into account (I201j) and the expression of Af (u) in (IllOp we get (I202[) . 

Again, a more precise notation for R m (u) and y? m (r c ) in theorem 3 would be 
R m ( r c,u) and y2 m (r c ,g), which emphasizes the explicit dependence on the critical value 
r c and the coupling constants g. 

8.4- The quartic model and the P}j equation 

Although the derivation of theorem 3 has been somewhat complicated, its applications 
are fairly simple. For example, we have already seen that the quartic model (|23l) at 
T c = g\j{A:g^) is a two-cut merging model of order m = 1 with r c = —giliA-g^). It also 
follows that if i = 2g2- Thus we have that the first coefficient u = di of a symmetric 
solution ( 1161!) of the corresponding string equation in the double scaling limit satisfies 
the P}j equation [22] 

g\u xx - 2g A (4g A u 3 + xu) = 0. (215) 
9. Concluding remarks 

In this paper we have presented a method to characterize and compute the large iV 
formal asymptotics of regular and critical Hermitian matrix models with general even 
potentials in the one-cut and two-cut cases. This method also leads to an explicit 
formulation, in terms of coupling constants and critical parameters, of the members of 
the Painleve I and Painleve II hierarchies associated with one-cut and two-cut critical 
models respectively. 
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As we pointed out in the introduction, the asymptotic form of the recurrence 
coefficient in the multi-cut case in general is not represented by (integer or fractional) 
power series in A -1 but involves a quasi-periodic dependence on N . Asymptotic 
expansions of the recurrence coefficients containing explicitly this quasi-periodic 
dependence on N have been considered by Bleher and Eynard [21] to prove the presence 
of the Painleve II hierarchy in a class of nonsymmetric models with merging of two cuts. 
It would be interesting to investigate if our method can be generalized to deal with 
this type of asymptotic expansions. This generalization would require new appropriate 
formulations of the continuum limits for both the recurrence coefficients r n ^ and s n ^, as 
well as for the generating function U n ^. A possible application could be the investigation 
of critical models featuring the birth of a cut [301 ED] , which do not seem to correspond 
to any conformal field theory and it is not clear if any integrable equation underlies the 
asymptotics of the associated recurrence coefficients in the double scaling limit. 
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